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Abstract 

Results on dynamical fluctuations of charged particles in the pseudorapidity 
space of central C-Cu interactions at 4.5 A GeV/c are performed in the transformed 
variables and using higher order scaled factorial moments modifyied to remove the 
bias of infinite statistics in the normalization. The intermittency behavior is found 
up to eighth order of the moments increasing with the order and leading to the 
pronounced multifractality. Two differed intermittent-like rises are obtained, one 
indicating an occurrence of the non-thermal phase transition, and no critical be- 
havior is found to be reached in another case. The observations may be treated 
to show different regimes of particle production during the cascade. Comparison 
with some conventional model approximations notes the multiparticle character of 
the fluctuations. The results presented can be effective in sense of sensitivity of 
intermittency to the hadronization phase. 
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1 Introduction 



At present intermittency [Q seems to be a well-founded fact observed in high-energy multiparticle pro- 
duction ^. This phenomenon, expressing as the power-law behavior, 

(F,)(xM'^% 0<^, <g-l, (1) 
of the g-order normalized scaled factorial moments (SFM), 

1 ^ (n^"^) 

m—l ^ ' 

("vertical" analysis identifies dynamical fluctuations. Here rim is the qth power factorial multi- 

nomial, Um {fim — 1) • • • {rim —9 + 1), with multiplicity rim in the mth bin from M ones, into which the 
space of the particles produced are divided. Average is taken over all events. 

The intermittency indexes, ipq, obtained via (|l]) are shown ^ ^ to be related to the fractal 
codimensions dq ||, |l^ , 

'Pq = {q-^)dq , (3) 

and correspondingly reflects the inner structure of fluctuations representing monofractal patterns with 
the unique dg, or multifractal ones with the hierarchy dq > dp ^ 1 > P- 

Applied to multihadron production processes, formation of such geometrical structures are pointed 
out ^1, |l2| to be a manifestation of one of two possible mechanisms. The first one, leading to the 
q-dependence of dq, is that one connects Q with self-similar cascade, rather a "non-thermal" (non- 
equilibrium) phase transition during the cascade than particle creation within one phase, e.g. hadronic 
[ p^ . The second scenario assumed at monofractality, is associated with thermal transition ^, e.g. 
from a quark-gluon plasma expected to be reached in central collisions of (ultra)relativistic nuclei [p^ . 

So study of intermittency decodes through fractality the geometrical and then thermodynamical 
features of high-energy multiparticle production. Note that the random cascading models being widely 
used as suitable to describe such processes (e.g., negative binomial distribution) JTst , lead in the most 
general form to the scaling (^ [Q , though they seem to be problematic to reproduce intermittency 
observed ^ |l8j. In general, the observations are still far from the qualitative explanation using 
existed particle-production codes, and origin of intermittency /fractality is a matter of debates 

Meanwhile, intermittency turns out to be very sensitive to the hadronic phase |^. Since the 
interpretation is mostly done in the approach based on the parton cascade models, it is rather difficult to 
explain the fact that partonic local fluctuations survive the hadronization process |l9| . In this context 
it is noticeable that "universality" of multihadron production in different type of reactions (from lepton- 
hadron to nucleus- nucleus) can be just a reflection of hadronization dynamics properties Thus, 
search for dynamical fluctuations and its possible treatment in "soft" process terms attracts considerable 
interest in reactions at intermediate energies. 

Important features of intermittency systematics behavior found in high-energy nuclear interactions 



[^-|26 are also already found to be manifested at relativistic energies. Besides purely hadronic model 
calculations showing intermittent structure due to the random cascade , enhancement of intermittency 
with beam energy decrease and weakening of the effect with increase complexity of reaction ^ 
are obtained. Recently, processes with low average multiplicity was discussed [ 3l[ |8| to dominate in the 
fractional SFM analysis. 



2 Experimental procedure 

The present paper deals with the study of the data obtained after processing the pictures of the 2m 
Streamer Chamber SKM-200 |Q equipped with a copper target. The chamber was installed in a 0.8 T 
magnetic field and it was exposed to the 4.5 A GeV/c ^^C beam at the JINR Synchrophasotron (Dubna). 
In data taking, the central collision trigger was used: the chamber was started if there were no charged 
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particles with momenta larger than 3 GeV/c in a forward cone of 2.4°. Details of the set-up design and 
data reduction procedure are described elsewhere |^ |3^. Systematic errors related to the trigger effects, 
low energy pion and proton detection, the admixture of electrons etc. have been considered in detail 
earlier Q and the total contribution does not exceed 3%. 

The scanning and handling of the film data were carried out on special scanning tables of the Lebedev 
Physical Institute (Moscow), using the method elaborated in ref. [3^ . The average measurement error 
in the momentum (Sp/p) was about 12%, and that in the polar angle measurements was (e^) ~ 2°. To 
search for dynamical fluctuations the charged particles in the pseudorapidity (77 = — lntan('i9/2)) region 
A77 = 0.2 — 3.0 (in the target rest-frame) were used, in which the angular measurement accuracy was 
not larger than 0.1 in the 77-units. The samples of 305 C-Cu events, which meet the above centrality 
criterion, have been selected with the average multiplicity of 27.2 ± 0.8 in the A77 under consideration. 

Earlier we have already analyzed fluctuations in the data presented and existence of non-statistical 
(dynamical) fluctuations was shown [^o| ^ |37| within intermittency approach and using multifractal 
analysis |^ . It is noteworthy that only applying the methods with statistical background suppression 
allowed dynamical nature of the fluctuations to be manifested, while the comparison between data and 
completely uncorrelated particle-production simulation fails to do that |]38| . It is significant in sense of 
the above discussed "hadronic origin" of the intermittency that analogous observation have been done at 
ultra-relativistic energies |Q (see also ref. ||39||). 

The study ^ also showed multifractality of the particle spectra along with multiparticle char- 
acter of the fluctuations being important to choice of the real particle-production model, as discussed. 
Multifractality have been also observed at ultra-relativistic heavy-ion collisions but obliged mainly to 
the two-particle correlations [|[ |j . Note that multiparticle contribution to the very-short-range correla- 
tions are directly observed in our data using the method of factorial cumulants |40[ as found in hadronic 
interactions recently ]4l|| . 

Before presenting the results, two important technical remarks should be made. First of all, in the 
previous investigations we have used "horizontal" analysis of the SFM taking into account the non-flat 
shape of the one-particle pseudorapidity distribution p{ri) to minimize the difference from "vertical" 
normalization (H) ^. Note that, though the "corrected" SFM were widely applied to analysis (in 
particular, at relatively small multiplicities) it implies the fluctuations to be bin-independent (the sum 
of the quotient in (0) transforms into the fraction of independent sums, X]m=i('^m)/X]m=i ^^^^ 



is, generally, a non-trivial assumption |4||5|]. In our data the slopes of the vertically-normalized moments 
are significantly greater the corrected ones at Srj < 0.5 [Q, while they coincides in collisions of ultra- 
relativistic ions p^ . 

To overcome the problem and, moreover, to compare the results observed in different experiments a 
new transformed variable, 



^ iv) = VIZl . , ' (4) 



have been introduced |4^, so that V is uniformly distributed in the [0;1] interval (^(^7) ~ const.). 
Due to the scale properties of the variables (^, the one-particle spectrum stretches in its central region 
(not significantly changing at the Ary-edges), eliminating losses from bin-splitting and, thus, allowing to 
observe higher-order moments. 

Another remark regards to the biased estimator of the SFM normalization, being the q-particle 
density function for uncorrelated production of particles in assumption of infinite statistics [0, 0, |^ . This 
sensitively influence the scaling law (|l|) for small bins. Note that a flattening of the moments for M > A/q 
is expected to be a reflection of the attainment of the correlation length Q . Really, measured SFM go to 
zero as the bin size aspires to the experimental resolution |^^, or much before because of the statistics 
limitations ("empty bin effect" Q ^). 

Recently, the modification of the method of the SFM have been proposed |Q to remove the bias in 
the normalization, and then the bias-free moments are defined to be , e.g. at "vertical" analysis. 



m— 1 
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where Nm is the number of particles in the mth bin in all M events, Nm = X]j=i("-ni)i- The property 
of the "vertical" and "horizontal" analysis to give the same results if the scaled variables are used 
(rim ~ {n)/M) seems to be also valid for the definition (||). 



3 Results and discussion 

In fig.l we show the log-log plots of the modified SFM (^) vs. number of bins of the space of "pseu- 
dorapidity" (^. Though by illustration reasons the dependencies are given only for some orders, viz. 
q — 3,5,6,8, they reflect the common peculiarities to be noticed. From the plots and the values of Lpq 
(see table) one can conclude that, besides considerable intermittent behavior, pronounced up to higher 
orders, there are the different power-like dependences in different intervals. Since at low orders {q = 2,3) 
this expresses as two sharply differed slopes: the values at M < 22 are about seven times less than ones 
at M > 23, this effect grows weak at g = 4 and 5. It is visible also that as the order increase an additional 
intermittent structure manifests at the range of small bins, but doesn't in fact survive when five and more 
particles are required to fill the bin and effect of small statistics becomes sizeable. 

However, regardless to the irregularities in the A/-dependence of the SFM, increase with the order 
up to higher moments. Moreover, in the M-interval of [7;17] and especially at 10 < Af < 17 strong 
intermittency is seen. 

It should be noted that the irregular behavior of the SFM was also observed in our previous investi- 
gation but using of the "ordinary" pseudorapidity variable strongly limited the order of the moment. 
Transforming 77-spectra into the uniform ones made this effect more pronounced and allowed to consider 
the SFM of high orders, where a new sub-structure is revealed. Note that strongly non-linear increas- 
ing of the log-log plots as presented have been observed in different interactions, from lepton-hadron to 
heavy ion collisions [Q, particularly at high-orders of the moments (e.g. [p2[) and/or at multidimensional 
analysis |2|, |2|, ||. 

As discussed, intermittency index, or rather codimension (^ dependence on g-order gives an infor- 
mation on the possible particle-production scenario. In fig. 2 we demonstrate dq as a function of q for 
different Af-intervals, within which a linear fit of the SFM plots (fig.l) is valid; the corresponding t^q- 
values are shown in the table. From the behavior of the dq obtained one can definitely conclude that very 
multifractal structure {dp > dq ,p > q) oi multiparticlc production is revealed, independent of the interval 
of M considered; no evidence for monofractality, and then for a second-order phase transition is seen. 
Let us to mention that decrease of the dq at high order is, in our opinion, connected with relatively small 
number of particles per event. Meantime, no considerable difference are obtained for some Af-intervals, 
namely for the couples of 7 < Af < 17, 10 < M < 22 and 4 < Af < 15, 2 < Af < 22. 

The multifractality, as noted above, lends support to find the condition for the non-thermal phase, not 
characterized by a thermodynamical behavior . As signal of the transition to this phase the existence 
of a minimum of the function 

Xq = -i 6 

q 

at a certain "critical" value of g = gc is required. In such a "spin-glass" phase the events are dominated 
by a few clots in projected distributions, while the "normal" phase is represented by the events with a 
bulk of peaks and holes |47, |4^. The minimum of eq. (H), as shown |47|, may be a manifestation of 



the fact that these two different phase mixed. Moreover, the phase transition can accompany (or occur 
inside) the branching process [|l^, ^. 

In fig. 3 the Xq are displayed as function of the (/-order for different Af-intervals with considerably 
distinguished dq (see fig. 2). A clear minimum at 7 < Af < 17 (10 < Af < 22) and 4 < Af < 15 
(2 < Af < 22) are seen, whilst the fits at 2 < Af < 28 do not exhibit such a feature. In our opinion, this 
is very intrinsic finding. 

Indeed, while sensitive behavior of the SFM in the region of 2 < Af < 22 are observing right up 
to eighth order (fig.l and the table) such a sharp difference of the dependence of Xq may be a signal 
of a non-trivial dynamical effectn. Referring to the non-thermal phase transition interpretation W9| 



The contribution of the statistics limitations is estimated to be small for these Af-intervals. 
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one should conclude that fits at different M-intervals lead to different cases of the intermittency: "weak" 
intermittency for monotone function (^) (large qc) and "stronger" intermittency when minimum is reached 
at g = 4 or 5. The stronger intermittency is related to the so-called "peak transition", i.e. transition 
at positive index qc, and the weak case contains the absence of the transition also. In our view, such 
an ambiguity is related to the cascade mechanism description that includes many steps irregardless to 
the change of the regime of particle production at different bin-averaging scales. This is, that apparently 
brings indetermination in treating of the intermittency observations e.g. in e+e^-annihilation | |49[ | and, 
meanwhile, being manifesting in our data is an extra proof of hadronic nature of intermittency observed 

. . . 

Noteworthily, only using of higher-order moments makes clear the hint to the non-thermal transition 
found earlier slightly visible up to fifth order |2^, |3^, |2^. Whether the minimum is taken place, the 
\ow-pt effect in the intermittency manifested at hadronic interactio ns Mlj leads to a clear minimum 
confirming responsible of hadronization for the effect of intermittency . The pt range for the stronger 
intermittency in our case was observed to be of 0.35-0.45 GeV/c like one have been already found 
searching for the maximum fluctuations It should be stressed also that since the minimum of the 
Xq ^ corresponds [|| to zeros of the fractal spectra, the system "frozen" at this point is no longer 
self-averaging and introducing of new observables is needed [^ . 

The shown importance of the cascade approach in the multiparticle production and the abovemen- 
tioned underlying of these processes to describe the hadronization attracts considerable interest to com- 
pare the predictions with the observations. Let us to limit oneself what have been earlier studied related 
to the SFM method, viz. negative binomial distribution (e.g. fl^ ), gaussian approximation Q and 
scale-invariant mass-splitting branching model . 

If the negative binomial distribution (NBD) is valid then the SFM are determined by the recurrence 
relation 

(^.) = (l + ^)(l + ^)---(l + 
where k is one of two NBD parameters should be independent of the Srj. The bin dependence of k 
expected from the intermittent rise at g = 2 (see table) is a reflection of the instability of the NBD 
ll). Besides, the NBD is a good fit of the data in the central 77-region only, and does not describe the 
"tails" of the multiplicity spectra, for which it transforms into the F-distribution being stable |Q . In the 
previous ref. we have shown that at the values of ipq closed to the experimentally obtained ones, the 
absolute values of the SFM approach each other when the shape of the ?7-spectra is accounted for. From 
this observation and the noted "central property" of the NBD the satisfactory agreement of the moments 
based on (0) with the experimentally measured ones in the "flat" spectra over the scaled variable ^ 
have been predicted. This is a fact well seen now from fig. 1 (x^ /NDF for the (Fq) are less than one 
standard deviation in all the A/-regions), may be excluding the case of g = 8 (e.g. x^/NDF « 6 at 
2 < M < 22). Meanwhile, the table shows that the NBD calculated slopes (pq yield the values markedly 
distinguished from the observed ones in the stronger intermittency case with possible non-thermal phase 
transition (7 < M < 17, 10 < M < 17), and trend to coincide each other at weak intermittency. In our 
opinion, this is because of input two-particle character of k calculated from (Q): flattening of the NBD 
prediction with q increase is visible at g > 6 (fig. 1). 

Such a difficult meat when the NBD is applied to the intermittency search makes it problematic to 
describe the data. Let us note that invalidity of the interpretation of multiparticle production at high 
energies in the NBD terms was also showed recently ^6, 17, 



The gaussian ( "log- normal" ) approximation (GA) predicts the relations for the slopes tpq to be defined 
as I 

^« = Y9('Z-1)- (8) 

In figs. 2 and 3 we show the GA predictions for the dq (||) and Xq (||) using the (/92-indexes at 
2 < M < 28 only. This is done by force of almost equal quantities of the d2 for all the cases of the 
Af -intervals considered. One can conclude the GA seems to be also hard to describe the intermittency 
effect especially when the phase transition is possible. In this sense it is interesting to note that the GA 
and the NBD approach both are found H] to describe the cascade with small number of steps, reflecting 
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invalidity of these approximations for the multiphcity asymptotics. The non-gaussian character of the 
correlations in multihadron production was also shown earlier ||, |l^ . 

The same result though less clearly is seen (fig. 2 and 3) if the simple scale-invariant cascade model 
proposed by Ochs and Wosiek is applied. In this case a modified power-law, 

(F.iSrj)) (X [gi6v)f^ , (9) 

was showed to be universal dependence for large class of models and of the available data of the one- 
and multidimensional SFM Here g{Srj) is the model-dependent function of {F2), so that 0^ = r^02. 
The conclusion agrees with the earlier observations p7[ | and shouts to one another obtained in hadronic 
interactions 

Thus all used model approximations, being based on the second order moments, indicate multiparticle 
character of the dynamical fluctuations at (possible) non-thermal (peak) phase transition. 



4 Conclusions 

The study of the intermittency phenomenon up to higher ranks of the scaled factorial moments of the 
pseudorapidity distributions of charged particles produced in central C-Cu interactions at 4.5 GeV/c 
per nucleon is performed. To eliminate the problem of the spectrum shape and to reach the higher 
multiple fluctuations the transformed variables are used. The unbiased modification is applied to the 
moments normalization to avoid a bias at small bins. The study shows existence of intermittent-like 
increase of the moments at all the orders considered, leading to a pronounced multifractality. The higher 
moment analysis allows to reveal two sharply differed increases with bin size, one of which indicates 
stronger intermittency and then non-thermal phase formation, and another one is usually connected 
with the weak intermittency when the critical g-order is not reached. While this comes from the fitting 
procedure at different bin intervals we are inclined to consider the fact obtained as a manifestation of 
different regimes of particle production at different scales of random cascading with a non-thermal phase 
transition inside. The results are compared to ones given in the assumption of a validity of the negative 
binomial distribution, gaussian approximation and scale-invariant mass-splitting branching model. The 
essential multiparticle character of the phase transition is indicated. Accounting possible hadronization 
influence on the dynamical correlations and search for the manifestations of a new matter formation 
the investigation of higher multiparticle fluctuations in nuclear collisions at intermediate energies gives 
evidence of further study of the effects observed. 
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Figure captions 

Fig. 1. The log-log plots of the modified scaled factorial moments (||) vs. the number of divisions. 
The curves present the NBD calculations (|^) and the straight lines show the least-squares fits at 10 < 
M < 17. 

Fig. 2. The codimensions dq versus g-order. 
Fig. 3. The Aq-functions (§). 



Table 1 



The intermittency indexes ipq compared to the NBD predictions. The errors present the covariance matrix 
estimators of the linear least-squares fits. 



q 


2 < M < 42 


2 < M < 28 


2 < A/ < 22 


4 < Af < 15 


2 


0.017±0.001 




O.OlliO.OOl 




0.009±0.001 




0.009±0.001 




3 


0.058±0.001 


0.059 


0.039±0.001 


0.035 


0.037±0.002 


0.029 


0.032±0.004 


0.030 


4 


0.130±0.002 


0.113 


0.087±0.003 


0.066 


0.092±0.004 


0.029 


0.079±0.009 


0.057 


5 


0.121±0.004 


0.180 


0.124±0.007 


0.106 


0.164±0.008 


0.088 


0.153±0.019 


0.092 


6 






0.113±0.015 


0.153 


0.233±0.013 


0.127 


0.234±0.035 


0.133 


7 










0.268±0.021 


0.172 


0.275±0.059 


0.180 


8 










0.207±0.031 


0.222 


0.207±0.091 


0.231 



q 


10 < A/ < 22 


7 < A/ < 17 


10 < Af < 17 


2 


0.004±0.002 




0.012±0.002 




-0.013±0.004 




3 


0.051±0.009 


0.020 


0.054±0.006 


0.035 


0.020±0.018 


0.00003 


4 


0.194±0.030 


0.039 


0.178±0.015 


0.069 


0.149±0.046 


0.0019 


5 


0.450±0.044 


0.062 
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